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CHAPTER  I 


INTRODUCTION 

Vibration  and  shock  have  become  important  fields  of  mechanical 
engineering  with  the  advent  of  high  speed  aircraft  and  missiles.  The 
equipment  or  "black  boxes"  installed  in  today's  air  vehicles  are  sub¬ 
jected  to  rigorous  vibration  inputs  and  sudden  shock  forces  and  must  be 
designed  to  operate  or  at  least  withstand  these  conditions.  Vibration 
can  be  defined  in  terms  of  displacement,  velocity,  and  acceleration 
while  shock  is  generally  defined  as  a  change  in  velocity  which  gives  rise 
to  high  accelerations.  Acceleration  is  the  value  most  often  desired 
because  the  destructive  force  is  proportional  to  the  mass  times  the 
indicated  acceleration.  Many  times  knowledge  of  the  destructive  force 
is  a  desirable  quantity  and  if  the  indicated  acceleration  cannot  be 
accurately  determined  the  calculated  force  will  be  of  little  or  no  value. 
Factors  such  as  mounting  techniques,  temperature  effects,  acoustical 
effects,  calibration  errors,  and  frequency  response  all  affect  the 
accuracy  of  the  recorded  environmental  data.  The  effect  of  the  frequency 
response  upon  the  accuracy  of  the  acceleration  data  will  be  the  main 
factor  discussed  in  this  paper. 

The  type  of  transducer  most  frequently  used  for  vibration  and  shock 
measurement  is  a  crystal  or  piezoelectric  accelerometer.  The  ability  of 
certain  crystalline  material  to  generate  a  charge  when  subjected  to  a 
physical  stress  has  been  known  for  many  years.  Some  of  the  comnon  crystal¬ 
line  materials  are:  Rochelle  Salt,  Quartz,  Barium  Titanate,  Lead 
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Zircon&te-Lead  Titanate  and  Lead  Metanlobate.  When  the  accelerometer  is 
subjected  to  an  input  a  mechanical  strain  is  applied  to  the  crystal.  The 
stress  in  the  crystal  causes  very  small  displacements  of  the  crystal 
domain  which  generates  a  charge  on  opposite  surfaces.  A  crystal  deflec¬ 
tion  of  three  micro -inches  may  produce  as  much  as  100  volts. 

The  input  to  be  considered  in  this  paper  will  be  a  shock  rather  than 
a  vibration  input.  The  shock  inputs  most  frequently  associated  with  the 
airframe  and  missile  field,  Insofar  as  environmental  testing,  are  a  10 
millisecond  square  wave  shock  pulse,  10  millisecond  half  sine  wave  shock 
pulse,  and  the  seven  millisecond  triangular  shock  pulse. 

The  ability  of  the  transducer  to  faithfully  reproduce  an  input  is 
related  to  the  frequency  response  of  the  accelerometer.  A  typical 
response  curve  for  a  crystal  accelerometer  is  shown  in  Figure  1. 


Fig.  1 

Frequency  Response  Curve  of  a 
Conventional  Crystal  Accelerometer 
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The  response  of  a  crystal  accelerometer  (Fig.  l)  with  respect  to  frequency 
has  a  useful  region  called  the  mid -band  range,  where  the  gain  is  independ¬ 
ent  of  frequency.  Above  and  below  this  mid-band  region  the  gain  is 
dependent  upon  the  frequency.  For  this  reason  it  is  always  desirable  to 
operate  the  accelerometer  at  frequencies  within  this  mid-band  range.  The 
question  then  arises  as  to  how  satisfactory  the  accelerometer  data  would 
be  for  certain  types  of  shock  excitation  where  the  frequency  content  of 
the  pulse  could  be  below  the  mid-band  range.  For  example,  if  the  frequency 
content  is  below  the  mid-band  range  (gain  leBS  than  unity)  and  the  input 
is  a  square  wave  the  output  of  the  accelerometer  will  not  be  a  "square" 
square  wave  (Fig.  2).  In  this  example  the  accelerometer  would  not  be 
faithfully  reproducing  the  input. 


Fig.  2 

Input  and  Output  Wave  Forms  for  Accelerometer  with  Frequency 
Response  lesB  than  Unity  in  the  Low  Frequency  Range. 

This  paper  will  present  a  method  of  calculating  the  error  of  a 
crystal  accelerometer  due  to  its  frequency  response  being  less  than  unity 
in  the  low  frequency  range.  The  inputs  considered  will  be  the  three  types 
of  shock  pulses  previously  mentioned. 


CHAPTER  II 


SHOCK  INPUTS  TRANSFORMED 
TO  rOSE  FREQUENCY  DOMAIN 


The  shock  input  is  a  function  of  time  (nonperiodic)  "F(t)M  and 
must  be  transformed  into  a  function  of  frequency  nF(f)",  This  is  neces¬ 
sary  to  determine  vhat  effect  the  frequency  response  of  the  accelerometer 
will  have  on  the  input  signal.  The  mathematical  tool  which  is  used  to 
transform  from  the  time  domain  to  the  frequency  domain  is  the  Fourier 
Integral  (Equation  l.l)1.  This  equation  is  analogous  to  the  Fourier 
Series  which  is  used  for  a  periodic  input  function. 


F(UJ)  „ 


P(t)  •  e 


-jUt 


dt 


£q.  1.1 


This  mathematical  tool  will  enable  the  engineer  to  determine  the  frequency 
content  of  any  nonperiodic  pulse,  by  this  it  is  meant,  the  number  of 
frequencies  which  must  be  summed  to  duplicate  the  pulse. 

Ihe  frequency  range  of  most  crystal  accelerometers  is  3-50,000  cps 
however  many  asgjlifiers  associated  with  the  accelerometer  have  a  range  of 
0-10,000  cyps.  It  is  recognized  that  there  are  a  number  of  suitable 

aiqpllfiers  whose  range  may  well  exceed  50,000  cps,  however  shock  measure¬ 
ments  are  primarily  associated  with  the  low  frequencies  and  this  paper 
will  deal  only  with  frequencies  between  0-10,000  cps. 


*V«n  Valkenburge,  M.  E.,  HETWORK  ANALYSIS .  Prentice-Hall  Inc., 
Englewood  Cliffs,  N.  J. ,  1955,  p.185. 
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Pig.  3 

Square  Wave  Shock  Pulse 


The  tine  function  of  a  square  wave  (Pig.  3)  is  given  by 

F(t)  -  Aq  -  |  <  t  <  |  E<1.  2.1 

and 

P(t)  =  0  everywhere  else 

where  "T"  is  the  period  of  duration  of  the  pulse.  Substituting  P(t)  -  AQ 
into  Equation  1.1,  P(f)  is  given  by 


T  T 
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2  K 

F(OJ)  - 

0 

~uT 

1  e  -  e 

i  i 

2  A„ 

F(co)  - 

0 

uy 

sin  ^2 

Eq.  2.2 

Substituting 

u  - 

2itf  and  rearranging  terms,  the  final  form  of  Eq. 

2.2  Is 

F(f)  -  A  E(l-  2.3 

0  2*77“  *  n 

where  ”f"  is  the  frequency  in  cycles  per  second  (cps). 

Equation  2.3  is  the  equation  for  the  frequency  spectrum  curve  for  a 
square  wave  shock  pulse  of  duration  "T".  The  graph  of  Equation  2.3  is 
shown  in  Fig.  A-l. 

HALF  SINE  WAVE  SHOCK  PULSE 


Fig.  4 

Half  Sine  Shock  Pulse 


Ibe  time  function  F(t)  for  a  sine  wave  (Fig.  4)  is  AQ  Sin 
between  0  and  T  and  zero  (o)  everywhere  else.  The  frequency  spectrum  for 
a  half  sine  wave  is  obtained  by  again  using  Equation  1.1  as  follows: 


F(  to)  - 


F(cj)  - 


F(to )  - 


F(cu)  - 


F(tu)  - 


F(to )  * 


F(u/)  - 


F(co  )  - 


F(lu)  - 


22 

2J 


f(|+W)(eJ" 


•  e-JUlT-l)  ♦  (f  -UJKe-J”  •  e-JtJT-l) 
J(f  -w)(f +W) 


An  t?  (Coa  x)  e"^  +  2U/j(Sin  x)  e'jgT  -  Ir 


-2  T 
2J 


J  ((*)  -W«) 


a  /  2x  -JtfT  2x\ 
-2  |  T  e  “  T 

2J  l  J  (|)"-US  / 
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F(W) 


F(u>) 


* 

_2_ 

T 


l  +  e 


-JWT 


A0  1  ♦  Coi  WT  .  ,1  sin 
T  **  - 


A  *  T 

F(w)  -  «*  -  (tVP  (1  +  Cos  WT  ■  J  Sin  WT) 

The  phase  shift  will  be  neglected  because  the  magnitude  of  the  frequency 
content  Is  the  primary  factor  in  this  paper.  Therefore , 


F(co) 

F(co) 

F(co) 

F(oo) 


VT 

*2  -  (,T  UJ)2 


V(l  +  Cos  WT)*  +  (Sin  WT)* 


A0*T 


-  (tu7F 

Aq  *tV2 

**  -  (tujP 


V5  +  2  Cos  WT 
Vl  +  Cos  WT 


A0*T 
-  (TUJ)^ 


Cos 


UJT 

2 


Eq.  2oU 


Substituting  U  *  2*f  into  equation  2.4  and  rearranging  terms  the 
final  font  is 

F(f)  -  i^TTZfT'' r»  *4*  2,5 

Equation  2.5  is  the  equation  for  the  frequency  spectrum  curve  for  a  half 
sine  vave  shock  pulse  of  duration  "T".  The  graph  of  Equation  2.5  is  shown 
in  Fig*  A-2* 
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TRIANGULAR  WAVE  SHOCK  PULSE 


Fig.  5 

Triangular  Shook  Pulse 


The  tine  function  of  a  triangular  wave  shock  pulse  (Fig.  6)  is 


F(t)  - 


0  i  t  <  T 


Eq.  2.6 


and 

F(t)  «■  0  everywhere  else. 

Substituting  Equation  2.6  into  the  transform  Equation  1.1  the 
resulting  frequency  spectrum  equation  is  given  by 


F(lu) 

F(uj) 

F(w) 

F(uo) 


J  v 


-JWT 
e  °  dt 


Ao  / 


-}U 


r  / 


ts  l  o'1* UI  ,  (_L.)S  (e-JWT  . 


of  1 4- *>- t  .1“»1 
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f(^i)  *  U'CT  +  J  T)  e’JWT  -  ‘O’] 

F(^)  ■  iJt  [(l7  +  -J  T)(Co8  UJT  -  J  Sln  UT)  • 

F(u j)  •  jj§  {$r  Kcos  u’1  -1)  +  T  Sln  UTl  +  J(T  Cos  WT  -  ■£  Sln  WT)} 


Again  neglecting  phase  shift 


F(u>)  = 

-2  /If 

wt  lTv 

F(u,)  - 

Ao  /  1 
wt  v  u T5 

F(^)  * 

a0^  r 

TJ5-f2  J  ■ 

F  (u.)  = 

y* ; 

W2  T2  ‘ 

Wz  T2 


-  Cos  UT  -  WT  Sin  WT 


(Cos  WT  +  WT  Sin  WT)  Eq.  2.7 


Substituting  ix)  =  2itF  the  final  form  is 

<c“ 


F(f) 


2«fT  +  2sfT  Sin  2*fT) 


Eq.  2.8 


Equation  2.8  is  the  equation  for  the  frequency  spectrum  curve  for  a 
triangular  wave  shock  pulse  of  duration  "T".  The  graph  of  Equation  2.8 
is  shown  in  Figures  A-3  and  A-4. 


CHAPTER  III 


INDUCED  ERROR  DUE  TO  THE  FREQUENCY 
RESPONSE  BEING  LESS  THAN  UNITY 

In  Chapter  II  the  three  types  of  shock  inputs  were  transformed  into 
the  frequency  domain  and  the  resulting  frequency  equations  were  plotted 
in  Figures  A-l,  A-2,  A-3,  and  A-U.  The  absolute  areas  hounded  by  these 
curves  are  used  to  determine  the  induced  error.  At  the  low  frequencies 
the  response  of  the  accelerometer  is  less  than  unity,  therefore,  those 
frequencies  in  the  input  signal  in  the  low  range  will  be  distorted  by 
the  transducer.  Ihis  means  the  absolute  area  under  the  frequency  spec¬ 
trum  curve  at  the  transducer  output  is  less  than  at  the  input.  The 
curve  at  the  output  is  the  product  of  the  frequency  spectrum  of  the  input 
pulse  times  the  frequency  response  of  the  accelerometer  (Fig.  7). 


Fig.  6 

Product  (Curve  C)  of  the  Frequency  Spectrum  (Curve  A)  for  a 
10  Millisecond  Square  Pulse  and  Frequency  Response 
(Curve  B)  of  the  Accelerometer 
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The  calculated  percentage  error  as  used  in  the  paper  ia  given  by 

$  error  ■  100 

The  transducer  that  vaa  chosen  as  a  typical  piezoelectric  accelero¬ 
meter  vas  the  Endevco  model  221?  manufactured  by  the  Kndevco  Corporation, 
l6l  East  California  Boulevard,  Pasadena,  California.  The  frequency 
response  (Fig.  8)  for  the  unit  is  based  on  a  total  capacitance  of  500 
mnfd  (accelerometer  and  cable)  and  an  amplifier  Input  liapedance  of  100 
megohms.  The  response  Is  flat  (unity  gain)  above  40  cps. 

For  all  curves  the  bounded  area  is  found  by  integrating  the  equation 
for  the  frequency  spectrum  curve  (2.3)  betveen  zero  and  10,000  cps.  This 
Integration  vas  accoiqpliehed  on  the  Bendix  0-150  digital  conputer  by 
using  Trapezoidal  Buie  vhich  is  the  numerical  integration  or  addition  of 
small  incremental  areas  along  the  frequency  axis.  The  general  form  is 

( -1-  y^-1)  (  a  f) 

where  and  y^+1  are  the  heights  above  the  frequency  axis  and  A  f 
is  the  frequency  increment  along  the  axis  (Fig.  7^. 

a 

Yu 


Fig.  7 

Frequency  Curve  Divided  into  Small  Incremental  Areas 
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Fig  8 

Frequency  ResponseaCurva  for  Endevco  Model  2213  Accalercwater 
(Capacitance  500ramfd  and  la pe dance  lOOmegohms) 
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ERROR  CALCULATION  RESULTS 

For  a  square  pulse  the  coeluted  area  is  1.501860  (Aq)  using  incre¬ 
ments  of  .25  cps  between  0.25  and  25  cps,  5  cps  between  25  and  3,000  cps 
and  20  cps  between  3,000  and  10,000  cps.  The  loss  in  area  due  to  the 
frequency  response  being  less  than  unity  below  25  cps  is  0.027238  (A^). 
This  is  the  difference  between  the  total  area  and  the  total  area  as  a 
result  of  distortion.  The  area  with  distortion  is  the  summation  of  the 
incremental  areas  times  the  gain  factor  for  each  respective  frequency. 
The  percentage  error  for  the  10  millisecond  square  wave  shock  pulse  is 

1 0.027238  (A  )  (loss  in  area)  \ 

100  1  1.501086  U0)  (total  area)  J  *  1>81^  * 

The  percentage  error  for  the  half  sine  wave  shock  pulse  of  10  milli 
seconds  duration  (using  the  same  approach  as  for  the  square  wave)  is 

0.03^355(A0)  (loss  in  area) 

0. 92572 (Aq)  (total  area) 


-  3.7m* . 


Likewise  the  percentage  error  for  a  triangular  wave  shock  pulse  of 
seven  milliseconds  duration  is 


^0.08t2888lU(AQ)  (loss  in  area)^ 
100  (  0. &3149729C A0)  (total  area)  . 


9.8988* 


CHAPTER  IV 


CONCLUSIONS 

The  accelerometer  error  due  to  the  frequency  response  of  the  Endevco 
unit  is 

1.8145J&  for  a  10  millisecond  square  wave  shock  pulse, 

3.7H1^  for  a  10  millisecond  half  sine  wave  shock  pulse, 
and 

9.8988^  for  a  seven  millisecond  triangular  wave  shock  pulse. 

The  error  for  the  square  wave  is  small  because  most  of  the  area  under 
the  curve  is  about  25  ops  where  the  gain  is  unity.  The  area  under  the 
curve  for  the  half  sine  wave  pulse  rapidly  decreases  for  increasing 
frequency  due  to  the  k^t2  term  in  the  denominator  of  Equation  2.5; 
therefore,  the  total  area  for  the  half  sine  wave  is  smaller  than  for  the 
square  wave.  However,  the  areas  below  25  cps  are  nearly  the  same  thereby 
making  a  larger  error  for  the  half  sine  wave. 

Bie  large  error  for  the  triangular  wave  is  due  to  the  fact  that  most 
of  the  area  in  the  low  frequency  range  where  the  gain  is  less  than  unity. 

While  the  results  of  this  paper  cannot  be  used  in  a  quantitative 
sense,  the  purpose  has  been  to  esghasize  this  particular  source  of  error. 
It  is  generally  neglected  in  practice  but  as  has  been  shown,  it  could 
represent  a  sizeable  amount  of  distortion. 
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Fig  A-l 

requency  Spectrum  for  a  Square  Ware  Shoe!  Pulse 
for  a  10  Millisecond  Duration 
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Fig  A-l 

Frequency  Spectrum  for  a  Square  V 
for  a  10  Millisecond  Dui 
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Fig  A-2 

Frequency  Spectrum  for  a  Half  Sine  Wave  Shock  Pulse 
for  a  10  Millisecond  Duration 
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Fig  A-3 

Frequency  Spectrum  for  a  Triangular  Wave  Sheet  Pulse 
for  a  7  Millisecond  Duration  between  5  and  2,200  ops 
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Fig  A-J* 

Frequency  Spectrum  for  a  Triangular  Wrve  Shook  Pulse 
for  a  7  Millisecond  Duration  between  0.25  and  6  cps 


